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Abstract. Mid infrared spectra provide a powerful probe 
of the conditions in dusty galactic nuclei. They variously 
contain emission features associated with star forming re- 
gions and absorptions by circumnuclear silicate dust plus 
ices in cold molecular cloud material. Here we report the 
detection of 6-8 /mi water ice absorption in 18 galaxies 
observed by ISO. While the mid-IR spectra of some of 
these galaxies show a strong resemblance to the heavily 
absorbed spectrum of NGC4418, other galaxies in this 
sample also show weak to strong PAH emission. The 18 
C\| | ice galaxies are part of a sample of 103 galaxies with good 
CD ■ S/N mid-IR ISO spectra. Based on our sample we find 
that ice is present in most of the ULIRGs, whereas it is 
weak or absent in the large majority of Seyfert s and star- 
\ burst galaxies. This result is consistent with the presence 
of larger quantities of molecular material in ULIRGs as 
opposed to other galaxy types. 

Like NGC4418, several of our ice galaxy spectra show a 
I maximum near 8 (ixa that is not or only partly due to 
■ PAH emission. While this affects only a small part of the 
Cn \ galaxy population studied by ISO, it stresses the need for 
^sO ■ high S/N data and refined diagnostic methods, to prop- 
erly discriminate spectra dominated by PAH emission and 
spectra dominated by heavy obscuration. 
| The spectral variation from PAH emission to absorbed 
. continuum emission near 8 /^m shows strong similarities 
"^h \ with Galactic star forming clouds. This leads us to be- 
O . lieve that our classification of ice galaxy spectra might 
reflect an evolutionary sequence from strongly obscured 
beginnings of star formation (and AGN activity) to a less 
' enshrouded stage of advanced star formation (and AGN 
C3 • activity), as the PAHs get stronger and the broad 8//m 
^ 1 feature weakens. 
• i— i . 
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1. Introduction 

Prior to the ISO mission, the mid-IR spectra of galaxies 
could only be studied from the ground in certain windows 
of reduced atmospheric absorption. For the far brighter 
Galactic sources less limitations applied, as these could 
be studied with airborne telescopes, like the Kuiper Air- 
borne Observatory (KAO: 1974-1995). Hence, before the 
advent of ISO, most of the prominent mid-IR ISM fea- 
tures had already been studied in some detail for Galactic 
sources, but not yet for extragalactic sources. Equipped 
with three mid-IR spectro(photo)meters, ISO has since 
enhanced considerably our knowledge of the mid-IR spec- 
tral properties of normal, starburst, Seyfert and Ultra- 
luminous Infrared Galaxies (ULIRGs). 

The mid-IR spectra of most galaxies are dominated by 
ISM emission features, the most prominent of which are 
the well-known PAH emission bands at 6.2, 7.7, 8.6, 11.3 
and 12.7 /im and atomic emission lines. The PAH features 
and the emission lines have been used qualitatively and 
quantitatively as diagnostics for the ultimate physical pro- 
ces ses po wering galactic nu clei (G enzel et al. 1998]; L utz et 
al. |1998|; Ri gopoulou et al. |1999| ; Clavel et al. |2000[ Helou 
et al. 200C; Tran et al. 2001 ). A broad absorption band 
due to the Si-0 stretching mode in amorphous silicates, 
centered at 9.7 /*m, is also commonly detected in galaxies. 
Since the center of the silicate absorption coincides with a 
gap between the 6.2-8.6/im and 11.3-12.8/im PAH com- 
plexes, it is not readily apparent whether a 9.7 /im flux 
minimum should be interpreted as the "trough" between 
PAH emission features or as strong silicate absorption, or 
as a combination of the two. 

In spectra observed towards heavily extincted Galactic 
lines of sight, like deeply embedded massive protostars, 
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a strong silicate feature is often accompanied by mid-IR 
ice absorption features due to molecules, frozen in grain 
mantles (Whittet et al. 1996 ; Keane et al. 2001 ). Among 
the simple ice molecules that have been identified in these 
grains are: H 2 0, C0 2 , CH 3 OH, CO and CH 4 ice. The 
first detections of interstellar ices in external galaxies were 
reported by Spoon et al. fl2000D and Sturm et al. fl2000D 
in the nuclear spectra of NGC4945, M82 & NGC253. At 
the low resolution of the ISO-PHT-S spectrophotometer a 
deep and broad 3.0 /im water ice feature, with the red wing 
extending far beyond the 3.3 PAH feature, is clearly 
visible in the spectrum of NGC4945. Also identified in 
the spectrum of NGC 4945 are absorptions due to 4.26 /im 
CO2 ice and the unresolved blend of 4.62 /im 'XCN' with 
4.67 /im CO ice and 4.6-4.8 /xm CO gas phase lines. 

The first detections of 6-8 /im ices in extragalactic 
sources were made in the strongly absorbed mid-IR spec- 



tra of NGC 4418 (Spoo n et al. |200l| ) and IRAS 00183- 
7111 (Tran et al. 2001). Especially rich is the spectrum 



of NGC 4418, which displays absorption features due to 
6.0 /im water ice, 6.85 & 7.3 yum hydrogenated amorphous 
carbons (HAC) and 7.67 /im CH4 ice, accompanied by a 
very deep 9.7 //m silicate feature. 

Encouraged by these findings we have searched our 
database of ^250 galaxies observed spectroscopically by 
ISO for galaxies showing similar 6-8 /im absorptions. In 
this paper we present the outcome of this search, which 
resulted in a sample of 18 galaxies showing evidence for 
the presence of 6.0 /im water ice. Sect. 2 describes the sam- 
ple. In Sect. 3, we discuss the complex interplay of 6.0 /im 
water ice and 5.25, 5.7, 6.25 /im PAH emission as well as 
the effect of rcdshift on the detectability of the blue wing 
of the water ice feature in ISO data. Sect. 4 describes the 
classification of the ice galaxies. In Sect. 5, we present the 
absorption and emission profile analysis. The results are 
discussed in Sect. 6. Conclusions are stated in Sect. 7. 

2. Observations 
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Fig. 1. A sequence of plots illustrating the effect of rebinning a pseudo high resolution spectrum (left panels) to the 
lower resolution of ISO-PHT-S and ISO-CAM-CVF (right panels). A redshift of cz=12000 is adopted and the grey 
scale represents the ISO-CAM-CVF rebinned spectra. The pseudo spectrum is composed of 6.0 /mi H 2 ice absorption 
and PAH emission bands found in the 5-7 fim range, superimposed on a flat continuum as shown in panel (a). In the 
remaining panels, (b, c, and d), the relative contributions of the H 2 absorption and PAH emission profiles are varied 
to investigate the behaviour of the spectral signatures. Note that we have discarded the 7.7 /im PAH feature in our 
toy model. This strong feature would start contributing to the spectrum at around 7 /im. 



2.1. ISO galaxy sample 

In order to study the mid-infrared spectral properties of 
galaxies we have created a database of some 250 ISO 
galaxy spectra. The database has been compiled from ob- 
servations performed with the PHT-S, CAM-CVF and 
SWS-01 instrument modes of ISO. The database com- 
prises: 



Nearly all ISO-PHT-S galaxy observations from a vari- 
ety of programs, except for galaxies at RA. 13-24 hrs 
from the "Normal Galaxy" program (P.I. G. Hclou). 
Four ISO-SWS-01 spectra of nearby galactic nuclei 
(Sturm et al. pOOOj ). 

ISO-CAM- CVF s pectra from the ZZULIRG proposal 
(Tran et al. 2001 ), as well as the a ctive galaxy spectra 
discussed by Laurent et al. ( 2000|) . 



By nature, this sample is not complete in any strict 
statistical sense. 



2.2. Selected sample 

We have defined a subsample of galaxies with good S/N 
spectra and sufficient wavelength coverage blueward of 
6.0/im. The criteria we applied were (S/N)6.5 Mm >8 and, 
for ISO-PHT-S in addition, cz>3000km/s. The resulting 
sample of 103 galaxies misses many of the observations of 
fainter ULIRGs and nearby galaxies in the ISO archive. It 
is dominated by Seyferts, which constitute nearly half the 
sample. 

Within this sample we found 18 galaxies show- 
ing 6.0 /im H 2 ice a bsorp tion (Figs.[| — 1|), including 
NGC 4418 (Spoon et al. [200l|) . We also looked for galaxies 
showing the opposite: no sign of water ice up to t icc <0.1- 
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Fig. 2. A collection of spectra showing strong 5.25 and 5.7/im PAH features (Class 6 galaxies; see Sect. 4). The four 
galaxies shown in the middle and right panels were observed using ISO-CAM-CVF, at low spectral resolution. The 
left panels show the spectrum of the Orion Bar, as observed at high spectral resolution by ISO-SWS (upper left) and 
degraded to ISO-CAM-CVF resolution (lower left). Note that the four galaxies shown represent the strongest cases of 
5.25 /im and 5.7/im PAH emission. The average galaxy may have far lower 5.25 /im/6.2 /im and 5.7/im/6.2/im PAH 
ratios. 



0.3 along the line of sight. Within our sample of 103 good 
S /N galaxy spectra, 28 galaxies fulfil this criterion. Exam- 
ples of ice-poor galaxies are shown in Figs.§, ||&[l0[ For 
the remaining 57 out of 103 galaxies, the quality of the 
spectra is either not good enough to detect the presence 
of a 6.0 /im water ice feature (a higher S/N criterion would 
have biased the sample against PAH-dominated spectra, 
which are weaker at 6-7 /im than continuum dominated 
AGN spectra), or to establish an upper limit better than 
0.3 for the amount of water ice. 

In the following Sections we show spectra for 32 
galaxies. These spectra were obtained using the ISO- 
CAM-CVF (Cesarsky et al. p.996[ ), ISO-PHT- S (Le mke et 
al. |1996D and ISO-SWS (De Graauw et al. |1996|) spec- 
trophotometers aboard ISO (Kessler et al. |l996j ). 

The ISO-CAM-CVF spectra of galaxies 100183- 
7111, 100188-0856 and Arp220 have been taken 
from Tran et al. ( pOOl] ). The ISO-CAM-CVF spec- 
tra of NGC253, NGC613, NGC 1068, M82, NGC 1365, 
NGC1808, NGC 4945, M83, Circinus and 123128-5919 



been reduced using standard routines of PIA[| versions 
8.1-9.0.1. The ISO-PHT-S spectra of NGC 4418 and 



have been taken from Laurent et al. (200C). 

The ISO-PHT-S spectra of galaxies Mrk334, NGC 23, 
101173+1405, NGC 828, 104385-0828, 105189-2524, 
MGC8-11-11, 106035-7102, UGC5101, Mrk231, 
Mrk273, Mrk279, 115250+3609, Arp220, 117208-0014, 
PKS 2048-57, 120551-4250 and 123128-5919 have 



NGC 4945 ha ve be en taken from Spoon et al. (|200lD and 
Spoon et al. ( [2000 ), respectively. 

The ISO-SWS spectrum of t he nu cleus of NGC 1068 
has been taken from Sturm et al. (200C), that of the Orion 
Bar from Peeters et al. (2002), and, finally, that of the 
nucleus of Circinus from Moorwood et al. ( |1996| ) . 

3. Ice absorption and PAH emission in the 
5-7)im region 

Before discussing the presence and strength of ice features 
in the sources of our sample, we illustrate the interplay of 
emission and absorption using a simple toy model. This 
is necessary since the 6 /im ice feature is not always seen 
against a well-defined continuum. The 'PAH' emission fea- 
tures at 7.7, 6.2, 5.7, and 5.25 /im may be present as well, 
leading to possible confusion of minima between emis- 
sions with absorption. In addition, the wavelength cover- 
age of the ISO data limits the accessible range, in partic- 
ular for low redshift objects observed with ISO-PHOT-SL 
(starting at 5.84 /im). Our model produces high resolution 
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spectra (R = X/AX ~ 1000) of the 5-7 /im spectral region, 
comprising an H2O ice absorption feature and PAH emis- 
sion bands superimposed on a continuum. Here, we de- 
scribe the various components of the model, the conse- 
quences of rebinning to the lower ISO-PHT-S and ISO- 
CAM-CVF resolution (i?- 40-90) and also explore the 
effect that extinction may have on the resulting spectra. 
In addition, we produce model-spectra at different red- 
shifts in order to investigate the effect that this has on 
the spectral signatures, in particular the relationship be- 
tween redshift and wavelength coverage, which is particu- 
larly important for ISO-PHT-S spectra. 

3.1. Components of the model 

The components chosen for the model closely mimic ab- 
sorption (ice) and emisson (PAH) features seen toward 
a variety of Galactic objects that have been extensively 
studied at high resolution. The spectra are modelled as: 
Fcontinuum x exp[— t(H 2 ice)] + PAHs. In order to allow 
easy comparison of the model generated pseudo-spectra, a 
false flat continuum is adopted. In Sect. 5 the flat contin- 
uum is replaced by continuum choices fitting the individ- 
ual observed spectra. Our model ignores any contribution 
from the PAHs to the overall continuum. This so-called 
PAH continuum can be seen in for instance the spectrum 



of the reflection nebula NGC 7023 (Moutou et al. |1998| ) as 
the continuum extending below the PAH features. Since 
the PAH continuum only serves to dilute the other con- 
tinuum and its associated ice feature, we ignore it here. 

Studies of Galactic star-forming regions, in particu- 
lar high mass embedded protostars, with the short wave- 
length spectrometer (SWS) of ISO show strong absorption 
features centred around 3.0 /^m and 6.0 /^m attributed to 
amorphous H2O ice. The spectral characteristics of H2O 
ice (whether pure or in various mixture s) hav e been well 
stu died in the laborator y (Ha gen et al. 1983 Hudgins et 
al. 1993| ; Maldoni et al. 1998 ). The spectral changes that 
occur during warm-up from 10 K to 120 K are irreversible. 
In the case of the 6.0 /im H 2 band, there is little spec- 
tral difference during warm-up from 10 K to 50 K though 
the band begins to broaden slightly and around 80 K this 
effect is becomes discernible. This broadening reflects an 
annealing of the amorphous structure. Further thermal 
cycling between 120 K and 10 K produces reversible ef- 
fects due to thermal shinking and expansion of the ice 
lattice which causes more minor variations in the width 
of the bands, particularly at 6.0 /im. Beyond 120 K the ice 
lattice transforms to polycrystalline ice. This is particu- 
larly noticeable in the 3.0 /im H2O band as sharp struc- 
ture begins to appear. Though the bulk of the interstel- 
lar 3.0 and 6.0 ^m absorption bands have been attributed 
to amorphous H 2 ice, there is also evidence of addi- 
tional molecules contributing to this feature in a num- 
berof sources (Dartois & d'Hendecourt 2001; Keane et al. 
2001 ). As a conseqence of this, the astronomical 6.0 /im 



feature varies somewhat from source to source and there- 
fore is not well suited to our simple model. Given the 
lack of 3.0 /im H 2 observational data for our sample, we 
have elected to represent the 6.0 /im ice absorption profile 
by a laboratory spectrum of pure amorphous H 2 ice at 
12 K. This profile has been shown to satisfactorily repro- 
duce the 6.0 /im H 2 feature observed toward a number of 



Galactic young stellar objects (Keane et al. 2001). The left 
panel of Fig.[l]a shows the profile of the 12 K H 2 feature 
(databa se of the Sackler Laboratory for Astrophys ics in 
Leiden: frrttp : // www . strw . leidenuni v . nl/ ~lab/ ) . The 
most striking aspect of this profile is the broad long wave- 
length wing which extends at least as far as 8.0 //m. There 
is little or no variation in the H 2 profile for tempera- 
tures < 50 K. For temperatures greater than this the broad 
wing becomes deeper and extends well beyond 8.0 /im. The 
right panel of Fig. |l|a shows the H 2 spectrum after it 
has been rebinned to the lower resolution of ISO-PHT- 
S (histogram) and ISO-CAM-CVF (grey line), assuming a 
cz=12000kms~ 1 . The profile is of course very similar but 
for the ISO-PHT-S spectrum there is incomplete pixel cov- 
erage of the short wavelength wing when compared to the 
ISO-CAM-CVF spectrum. The dominant PAH emission 
features (6.2, 7.7, 8.6, 11.3 /im), observed toward a variety 
of Galactic sources, show a wide degree of variation in the 
detailed profile shapes and relative intensities (Peeters et 
al. 2002). On the other hand, observations of PAH emis- 
sion bands in external galaxies reveal very little variation 



in the profiles (Rigopoulou et al. 1999; Helou et al. 2000) 



Since one of the main aims of the model is to investi- 
gate the effect of rebinning spectral features to the lower 
resolutions of ISO-PHT-S and ISO-CAM-CVF, adopting 
extra-galactic PAH emission profiles is inappropriate. In- 
stead, high resolution Galactic PAH emission band pro- 
files are adopted. The ISO-SWS PAH emission spectrum 
of the Orion bar PAH is chosen for the model as this is 
a good represe ntatio n of a region of active star-formation 
(Peeters et al. fc002|) . The left panel of Fig.Qa shows the 
high resolution spectrum of the Orion Bar. Apart from 
the typical 6.2 /im PAH CC stretching mode, the observed 
spectrum also shows evidence for two weak PAH emission 
features near 5.25 /urn and 5.7 /im (left panel of Fig.|l|a). 
The nature of these weak PAH features is not very well 
studied but they are believed to be combinations or over- 
tone bands involving C-H bending vibrations (Allaman- 
dola et al. 198£). The strength of these weak bands rel- 



ative to the strong 6.2 /im band is variable. Of a sample 
of 35 compact HII regions studied by Peeters et al. (in 
prep.) roughly 10% of the sources show evidence for these 
weak PAH bands. Consequently, because of their variable 
strength, their presence is not always assured in spectra 
showing the dominant PAH bands and there is no rule of 
thumb defining an expected ratio for the intensity of these 
weak PAH bands to the 6.2 /im band. The Orion bar is 
among the Galactic sources with the strongest 5.25 and 
5.7 /im features relative to the 6.2 /im feature. Compar- 
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ison with high S/N CAM-CVF spectra of the brightest 
nearby starbursts suggests it to be fairly representative 
for those objects however (Fig.||, confirmed by SWS spec- 
troscopy of M82 and NGC253 (Sturm et al. |2000| )), We 
hence adopt the Orion Bar as a template for the 5.25, 5.7, 
and 6.2 /im PAH features. The PAH spectrum in Fig. [l] was 
obtained by subtracting a spline continuum from its spec- 
trum. Fig. |l| then suggests that in the presence of strong 
PAH emission, a shallow 6.0 /im ice absorption of less than 
about 10% of the peak height of the 6.2 /im PAH feature is 
very difficult to discriminate from the minimum between 
the 5.7 and 6.2 /im PAH features, especially for the lim- 
ited wavelength coverage of ISO-PHT-S. Finally, the 7.7 
and 8.6 /im PAH bands have been clipped from the spec- 
trum as they will be dealt with separately (Spoon et al. 
in prep.). 

The right panel of Fig. |l|a shows the Orion-Bar PAH 
features rebinned to the lower resolution of ISO-PHT- 
S (histogram) and ISO-CAM-CVF (grey line). In rebin- 
ning, the peak height of the 6.2 /im feature is reduced, 
more for ISO-CAM-CVF than for ISO-PHT-S. In the 
ISO-CAM-CVF spectrum, the weak PAH bands are still 
present, though they are flattened and broaded as com- 
pared to the high resolution spectrum. Only the peak of 
the 5.7/xm weak PAH feature is present in the ISO-PHT- 
S spectrum since at a cz of 1 2000 km s^ 1 the first pixel is 
at 5.6 /im. This is discussed further in Sect. 3.3. Finally, 
the location of the water (i.e., in front, behind or mixed 
with the PAHs) had little effect on the resulting shape of 
the spectral features, in particular on the profile of the 
6.2 /im feature. We hence ignored extinction on the PAHs 
in our toy model. Detailed fits to the extragalactic spec- 
tra are presented in Sect. 5, in which the flat continuum 
is replaced by continua determined individually for each 
galaxy. 

3. 2. Effects on the model of varying the component 
contributions and extinction 

To investigate the interplay of the H2O absorption feature 
and the 6.2 /im PAH emission bands, the relative con- 
tribution of each component is varied. This is shown in 
Fig.|l|b-d, where the left panel is the high resolution in- 
put spectrum and the right panel represents the rebinned 
low resolution pseudo spectra of the model. The dotted 
line indicates the flat continuum. Panel (b) illustrates the 
effect of a strong H2O band on the PAH emission bands. 
Regardless of the fact that the 6.2 /jm features lies within 
the long wavelength wing of the H2O feature, its profile 
is preserved (left panel). The 5.7 /xm weak PAH band, on 
the other hand, is lost within the short wavelength wing, 
though there is a slight hint of a peak at 5.7/zm. The 
5.25 /im PAH feature is still present as the H2O profile 
does not reach as far as this short wavelength position. 
Rebinning to the low resolution of ISO-PHT-S and ISO- 
CAM-CVF (right panel) shows that the H 2 feature, in 
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Fig. 3. Panel (a) shows a mid-IR extinction law of the 
Galactic centre derived from hydrogen recombinatin lines 
(Lutz 199S| ). Illustrated in panel (b) is the effect of ex- 
tinction on the model spectra for different amounts of A v , 
assuming the above extinction curve. The effect is most 
severe for the 8.6 /im band as all profile information be- 
comes lost within the 10 /im silicate feature. 



combination with the 6.2 /im PAH emission band, domi- 
nate the model spectrum and there is little if no evidence 
for the weak PAH bands in the ISO-PHT-S spectrum. 
Only the ISO-CAM-CVF spectrum reveals the presence 
of the 5.25 /im weak PAH band, due to the wider wave- 
length coverage. As the strength of the H2O absorption 
band is reduced (Panels (c) and (d)), the 6.2 /im feature 
starts to dominate the model spectra and now only weakly 
sits within the long wavelength wing of the H 2 feature. 
Also, the spectral signature of the 5.7 /im weak PAH band 
starts to become pronounced and in Panel (d) the struc- 
ture is clearly distinguished lying above the flat contin- 
uum. The crucial result of this study is that the profile of 
the 6.2 /im PAH band is little affected by the presence of 
an H2O ice absorption feature. In addition, it is evident 
from this analysis that the weak PAH bands do not influ- 
ence the overall characteristics of the structure seen in the 
model spectra. However, the weak PAH emission features 
may enhance/mimick the effects of weak ice absorption. 

The possible effect that extinction may have on the 
spectral features of the pseudo spectra were also investi- 
gated. In light of the fact that the Galactic Center has 
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been well studied and also because of the homogeneous 
nature of the foreground extinction, the extinction law 
derived by Lutz ( |1999| ) for the Galactic Centre is adopted 
(Fig.||a). We have also investigated the effect of other ex- 
tinction laws (Tran et al. |2001 ; Draine 1989] ) but they are 
substantially the same and only the effects of the Lutz 
Law are discussed in detail here. The extinction is applied 
as if the continuum+PAH+ice combination were behind 
a column of material. Though the situation is less clear 
for external galaxies, this approach allows for a qualita- 
tive assessment of the effect of extinction on the model 
spectra. In order to clearly demonstrate the effect of ex- 
tinction, the 7.7 and 8.6 /im PAH bands of the Orion Bar 
have been included. The top plot in Fig.J^b shows the un- 
extincted (A v =0) model spectrum. As the degree of ex- 
tinction increases from A v —Q to A v =25 the 8.6/xm PAH 
band is " eaten away" by the short wavelength wing of the 
10 /im silicate feature. Eventually, the depth of the silicate 
band is so great that the 8.6 /im feature is completely lost 
within this band. The same effect, but with less sever- 
ity, is noticeable for the 7.7 /im emission feature. In the 
case of the 6.2 /im band, the peak intensity decreases as 
the degree of extinction increases, but information on the 
profile shape remains intact and distinguishable. Thus, the 
spectral characteristics of features shortward of 7 /zm are 
unaffected by the effects of extinction. Furthermore, due 
to the sensitivity of the features beyond 7 /im to the degree 
of extinction they have been neglected in the model. 

3. 3. Profile appearance as a function of redshift 

We have also investigated the effects of redshift com- 
bined with the limited spectral coverage on the appear- 
ance of the spectrum. A redshift (cz) range of Okms -1 
to 20 000 km s" 1 was chosen as this is representative of 
the cz range of the observed galaxies (Table [l]) . Fig. |] il- 
lustrates the spectral variation as a function of redshift. 
The high resolution spectrum shown in panel (b) of Fig. [j] 
was adopted as the input to the rebinning routines and 
as before the histogram and the grey line represent the 
ISO-PHT-S and ISO-CAM-CVF rebinned data, respec- 
tively. Since the spectral coverage of ISO-CAM-CVF ex- 
tends down as far as 5 /im there is not much variation 
in the position of the first pixels as a function of red- 
shift. Hence in Fig. [| only one ISO-CAM-CVF spectrum is 
shown. For redshifts of 6000 km s -1 or greater, the 6.0 /im 
H2O absorption profile is readily seen in the model ISO- 
PHT-S spectra. On the other-hand, the presence of H 2 
is not so obvious in ISO-PHT-S spectra for redshifts less 
than 6000 km s _1 . As the redshift decreases the pixel cov- 
erage of the short wavelength wing of the H2O absorp- 
tion feature diminishes and only one or two pixels begin 
to rise at the shortest wavelengths. Consequently, there 
is not much real evidence for the blue wing of the H2O 
band even though the original high resolution model spec- 
trum clearly shows strong H2O ice absorption (panel (b) 
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Fig. 4. A comparision of the spectral appearance of panel 
(b) from Fig. [j] rebinned to ISO-PHT-S resolution as a 
function of redshift. At low redshifts the presence of 6.0 
/im H2O ice absorption is not very convincing. However, 
as the redshift is increased the starting position of the 
first pixel shifts to lower wavelengths resulting in more 
of the 6.0 /im H2O profile being clearly discernible. Since 
the wavelength coverage of ISO-CAM-CVF (grey line) is 
significantly larger than that of ISO-PHT-S, there is not 
much variation in wavelength of the first pixel with red- 
shift. 

in Fig. Q). In the case of the 6.2 /xm PAH emission band, 
the only notably change is that the rebinned structure of 
the profile peak varies with redshift. 

4. Classification 

Although likely the spectra of galaxies represent a contin- 
uous distribution with variable absorption and emission 
components, we have classified the galaxies in our sample 
into six categories, based on their mid-IR spectral proper- 
ties. The first three catagories comprise the galaxies with 
evidence for the presence of 6.0 /im water ice absorption 
in their spectra: 

— Class 1: These galaxies (Fig.||) exhibit a 6.0 /im 
H2O ice absorption feature, without contamination 
by 6.2 /im PAH emission. Galaxies of this type also 
show other 6-8 /im absorption features - most com- 
monly the 6.85 /im feature, attributed to HAC. Note 
that all galaxies in this class have a broad flux peak at 
~8/im. Judging from their complete 2-200 /im SEDs, 
all Class 1 galaxies have a deep 9.7/im silicate feature. 
For the two galaxies for which 2-5 /im ISO-PHT-S data 
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is available, NGC4418 and IRAS 15250+3609, there is 
a sharp increase in flux between 5 and 6 /im. 

— Class 2: These galaxies (Fig.^|) show 6.0 /im H2O ice 
absorption features, partially fillcd-in by weak 6.2 /im 
PAH emission. Other 6-8 /im absorption features are 
not always as apparent. Note that like in Class 1, all 
galaxies in this class have a broad 8 /im feature. 

— Class 3: These galaxies (Fig. exhibit a 6.0 /im 
H2O ice absorption feature, largely filled in by strong 
6.2 /im PAH emission. No other 6-8 /im absorption fea- 
tures have been found. All galaxies in this class show a 
7.7/im PAH emission peak, clearly narrower than the 
broad 8.0 /im feature seen in Class 1 & 2 galaxies. 

Table|l] lists the classification for all 18 ice galaxies. 
For completeness we also list galaxies with 3.0 /im ice 
absorption in their nuclear spectra: NGC253, M82 and 
NGC4945 (Sturm et al. pOOO} Spoon et al. |2000|). Only 



strong as Ti co =0.3. We note that among the Class 6 galax- 
ies, both NGC 253 and M82 do have a weak 3.0 /im water 



for NGC 4945 has the usually far weaker 6.0 /tm feature 
also been detected - at similar strength as the 3.0 /tm one. 
This apparent inconsistency may be due to aperture size 
effects (the 3.0 /im feature was measured by ISO-PHT-S; 
the 6.0 /im feature by ISO-CAM-CVF) or due to dilution 
by less obscured emission, filling in the water ice features. 

Within our sample of 103 galaxies with good S/N spec- 
tra and sufficient wavelength coverage we also recognize 
three other classes of spectra: 

— Class 4: These galaxies (Fig.|J) show a smooth fea- 
tureless continuum in the 5-6 /im range as well as defi- 
nite signs of 9.7 /im silicate absorption. The presence of 
silicate absorption is most easily determined by inter- 
polating the 6 /im ISO and the 12 /im IRAS flux points 
(or in case of ISO-CAM-CVF, using the 13-15 /im con- 
tinuum instead). 

— Class 5: These galaxies (Fig.|To|) show a smooth fea- 
tureless continuum in the 5-6 /im range without clear 
signs of 9.7 /tm silicate absorption. The absence of a 
clear silicate absorption feature is most easily deter- 
mined by interpolating the 6 /im ISO and the 12 /im 
IRAS flux points (or in case of ISO-CAM-CVF, using 
the 13-15 /tm continuum instead). 

— Class 6: These galaxies (Fig.|J) clearly show 5.25 and 
5.7 /tm PAH features, without clear signs of 9.7 /tm sili- 
cate absorption. The absence of a clear silicate absorp- 
tion feature is most easily determined by interpolating 
the 6 /tm ISO and the 12 /tm IRAS flux points (or in 
case of ISO-CAM-CVF, using the 13-15 /tm continuum 
instead) . 

Table U lists the observed mid-IR features for the 13 
out of 28 Class 4-6 galaxies displayed in this paper. It is 
hard to quantify consistently the limits on ice absorption 
in these three classes because of the presence of PAH fea- 
tures. For Class 4 and 5 we estimate Ti ce <0. 1-0.3. For 
Class 6, in all spectra the presence of weak 5.7 /tm and 
strong 6.2 /tm features may easily mask ice features as 



ice feature (Sturm et al. 2000), but the 6.0 /tm ice limit is 
T icc =0.3 for both galaxies. 

The total number of ISO galaxies classified into Classes 
1-6 is small: just 46 out of ~250 galaxies in our sam- 
ple. The combination of both a good coverage of the 5.5- 
6.0 /tm range and a good S/N over 5.5-6.5 /tm proves to 
be a hard requirement to meet. With its superior sensi- 
tivity and better coverage of the 6/tm region, SIRTF can 
be expected to refine our classification, find new members 
and provide a larger sample for statistical analyses. 



5. Absorption and emission profile analysis 

5.1. Model fits 

The model described in Sect. 3 is now applied to the Class 
1-3 sources presented in Sect. 4. The flat continuum of 
Sect. 3 is replaced by the 5-8 /im continua derived individ- 
ually for each source (thick grey lines in Figs.||-^). For 
Class 2 and 3 sources, the continua were determined by 
interpolating the ISO-PHT-S short wavelength spectrum 
(2-5 /tm) linearly or logaritmically (depending on the SED 
shape) to the long wavelength IRAS data, ignoring the 
presence of PAH emission features. For Class 1 sources we 
followed another recipe. In the absence of contamination 
by PAH emission bands, we assumed the 8 /tm flux peak to 
represent the local 8/zm continuum (see Sect. 6 for a dis- 
cussion on the nature of the 8 /tm flux peak). For the three 
sources observed with ISO-CAM-CVF we then derived the 
5-8 /tm continuum by interpolating to the observed 5 /tm 
flux. For NGC 4418 and 115250+3609 we interpolated to 
the reddest 2-5 /tm ISO-PHT-S data point, assuming a 
continuum shape similar to that of 100188-0856. Note 
that for all Class 1 sources we included a 6.85 /tm ab- 
sorption feature du e to h ydrogenated amorphous carbon 
(HAC; Furton et al. 1999 ). A (strong) 6.85 /tm absorption 
feature is commo nly ob served towar ds G alactic lines of 
sight (Chiar et al. |2000| ; Keane et al. fcOOlp , 

Fig. 11 displays the model fits to all the sources, ex- 
cept for the three sources discussed below (Fig. fL2|). The 
grey histogram represents the model and in all cases the 
model spectra match very well the observed data. For two 
sources, Arp 220 and NGC 4945, the ISO-CAM-CVF spec- 
tra have also been included (thin line). Reading from left 
to right, of the first 5 panels in Fig. [ll] a HAC absorption 
feature has been included for four of the sources. This 
profile matches very well the strong absorption feature 
seen near 6.85 /im in these sources. In the case of 1 00188- 
0856, absorption between 6.4 and 7.4 /tm is not consistent 
with the 6.85 /im absorption profile and therefore only the 
6.0 /tm H2O profile is modeled. We are not aware of any 
other valid or likely candidate that can give rise to the 
absorption observed in 1 00188-0856. For the remainder of 
the sources, this feature has been omitted from the model. 
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Fig. 5. Class-1 mid-IR galaxy spectra. Galaxies in Class 1 exhibit a clear 6.0 /im ice absorption feature, without clear 
evidence for the presence of a 6.25 /im PAH feature. All sources in this class have in common that they have a broad 
continuum flux peak near 7.7 fim, clearly broader than a 7.7 (im PAH feature. The panels on the left show the 5-12 /im 
spectra on a linear scale, while the panels on the right show the same spectra on logarithmic scale, with smoothed and 
rebinned 2-5/mi ISO-PHT-S spectra and 12-100 /im IRAS fluxes added. All fluxes have been k-corrected. The thick 
grey lines show what the 5.5-8.0 /im continuum would be like if no ice absorption were present. The vertical error 
bars at 5.8 and 11 /im denote the pixel to pixel flux error (la, so not Her), while the horizontal error bars denote the 
gradually changing ISO-CAM-CVF spectral resolution (R~40). 



A qualitative comparison of the panels of Fig. |TTj, reveals 
that the observed features at 6.0 /im and 6.2 /im are suc- 
cessfully reproduced for all sources. 

For three of our ice galaxies (Mrk231, NGC828, 
Mrk334) proving the presence of 6.0 /xm ice absorption - 
the very key to membership of Classes 1-3 - turned out to 
be non-trivial. For each of these galaxies (see Fig. |l2j) the 



number of pixels in the blue wing is small, as is the depth 
of the ice feature. In all cases continua exist which do not 
require the presence of 6.0 /im ice absorption. However, 
such continua, which are defined here as the superposi- 
tion of hot dust plus PAH continua, must give a good fit 
in the 6-8 /im range (i.e., any emission above or below the 
continuum level should be accounted for by known emis- 



10 



H.W.W. Spoon et al.: Ice features in the mid-IR spectra of galactic nuclei 




10.0 r 



1.0 r 



0.1 



10 



100.00 r 



10.00 



1.00 



0.10 r 



0.01 



H 1 1 1 — I I I I 



05189-2524 

H 1 1 1 1 — I I I I — 



H 1 1 1 — I I I I 



Mrk231 



H 1 1 1 1 — I I I I 





20551 -4250 



9 10 11 12 



10 



20 



60 100 



Rest wavelength [/irm] Rest wavelength [/im] 

Fig. 6. Class-2 mid-IR galaxy spectra. Galaxies in Class 2 exhibit a clear 6.0 /im ice absorption feature, with clear 
evidence for the presence of a 6.25 /im PAH feature. All sources in this class have in common that they have a broad 
continuum flux peak near 7.7 yum, clearly broader than a 7.7 (im PAH feature. The panels on the left show the 5-12 fim 
spectra on a linear scale, while the panels on the right show the same spectra on logarithmic scale, with smoothed and 
rebinncd 2-5/im ISO-PHT-S spectra and 12-100 /im IRAS fluxes added. All fluxes have been k-corrected. The thick 
grey lines show what the 5.5-8.0 /im continuum would be like if no ice absorption were present. For Arp220 both the 
ISO-PHT-S (black) and ISO-CAM-CVF (grey) spectrum are shown in the left panel. The vertical error bars at 5.8 
and 11 /im denote the pixel to pixel flux error (la, so not ±1<t). 



sion and absorption features) and the continuum should 
join smoothly to the long wavelength (>12/im) data. In 
the case of Mrk 231, there is a preference for the fit requir- 
ing the presence of water ice (Fig. 12, left panel). This fit 



adopts a much more realistic continuum (cf. Fig. 6) and 
it provides a more realistic value for the 6.2 /xm/7.7/im 
PAH ratio than the fit not requiring the presence of water 



ice (Fig. fl2[ right panel). Moreover, the fit without wa- 
ter ice would result in a very broad 7.7 /zm PAH feature 
which does not resemble PAH features in Galactic sources. 
However, the total amount of ice absorption is difficult to 
quantify. A slightly lower but still permissible continuum 
can reduce the ice optical depth by a factor 2. For the 
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Fig. 7. Class-3 mid-IR galaxy spectra. Galaxies in Class 3 exhibit a narrow 6.0 /im ice absorption feature, partially 
filled in by the wings of the 6.25 fim PAH feature. All sources in this class have in common that the flux peak near 
7.7 /zm has the typical PAH width. The panels on the left show the 5-12 /im spectra on a linear scale, while the panels 
on the right show the same spectra on logarithmic scale, with smoothed and rebinned 2-5/rai ISO-PHT-S spectra 
and 12-100 (im IRAS fluxes added. All fluxes have been k-corrected. The thick grey lines show what the 5.5-8.0 ^m 
continuum would be like if no ice absorption were present. For NGC 4945 both the ISO-PHT-S (black) and ISO-CAM- 
CVF (grey) spectrum are shown in the left panel. The vertical error bars at 5.8 and 11 /xm denote the pixel to pixel 
flux error (la, so not ±ler). 
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other two galaxies (NGC828, Mrk334) the fits with and 
without water ice are equally viable. 

5. 2. Derived physical parameters 

Table | summarizes the physical parameters of the de- 
tected features. The range in optical depths of the H 2 
absorption band is large, extending from small depths 
of T icc =0.13 (105189-2524) to very deep features with 
r icc =1.3 (UGC5101). The optical depths of the H 2 
band were converted to column densities by dividing the 
peak optical depth by the intrinsic peak strength (4.2 
xl0~ 20 cm 2 molecules -1 ; Hagen et al. 1983| ). The derived 
column densities are comparable to H2O column densi- 
ties determined for Galactic molecular clouds (Keane et 
al. |2001| ). All sources in Class 1, except for 100188-0856, 
also show evidence for additional absorption features cen- 
tred at 6.85 /im, 7.3 /im, 7.67 /im. The 6.85 /im feature re- 
sides in the long wavelength wing of the H2O band, and 
hence the depth of this feature is more accurately deter- 
mined by dividing out the continuum plus H2O feature 
rather than just the continuum. For Class 1 galaxies, ex- 
cept for 100183-7111, the depth of this feature is very 
strong and is greater than the corresponding H 2 optical 
depths. In Galactic dense molecular clouds, on the other 
hand, the water ice absorption is always the stronger of 
the two (Keane et al. 2001). The origin of the 6.85 /im 



feature is unknown. The presence of a 7.3 /im feature in 
some spectra suggests it is carried by HAC residing in the 
diffuse ISM , very much like the Galactic Centre (Chiar 
et al. 2000 ). An absorption feature at 7.3/zm is present 
in the sp ectra of 1 15250+3609 and NGC4418 (Spoon et 
al. 2001 ). On the other hand, many Galactic young stel- 
lar objects show a strong 6.85 /im feature due in part to 
an unidentified ice component (Keane et al. 2001 ). Ob- 
servations of the 3 /zm region could distinguish between 
these two possibilities because HAC material presents a 
strong 3.4 /im absorption feature which is absent in Galac- 
tic YSOs (Pendleton et al. [l994t Pendleton & Chiar |1997| ). 
A 3.4 /im absorption feature observed towards a few extra- 
galactic sources has been tentatively attributed to HAC 
(Wright et al. |1996| ). Finally, one source, NGC 4418, has an 
absorption feature at 7.67/im which has been attributed 
to methane ice (CH 4 ; Spoon et al. pOOl] ) . 

Since the 6.2 /im feature lies, in general, within the 
wing of the H 2 absorption band and hence slightly be- 
neath the continuum, it is necessary to subtract a combi- 
nation of the continuum plus H 2 from the observational 
data. After subtraction, the intensity (Wcm -2 ) of the 
6.2 /im band is determined between 6.0 /im and 6.55 /im. 
The derived 6.2 /im PAH intensities are listed in Table|l|. 
Note that for models in which the water ice is mixed with 
or in front of the PAH emitting region, the intrinsic 6.2 /im 
PAH intensities would be greater. 

An estimate of the 9.7/im silicate optical depth (t s h) 
for each of the sources is also given in Table 0. In order 



to be able to compute this quantity the 9.7/im continuum 
was interpolated from the 5-8 /im continuum and the long 
wavelength data (>12 /im). For some sources the peak ab- 
sorption of the silicate feature is saturated and the true 
depth of the band is thus uncertain. In these cases fitting 
the wings of the feature can be attempted. Given the pos- 
sibility of additional sources of emission along the line of 
sight, however, we prefer here to just state lower limits for 
r 8 u. The column density of hydrogen is computed from the 
9.7/im silicate depth by assuming a Galactic conversion 
factor:- N H = ^- x ^ x t sU = r sU x 3.5 x 10 22 cm -2 ; 



(Roche & Aitken |198J, Bohlin et al. |1978[ ) 



6. Discussion 

Within our sample of 103 galaxies with good S/N and suf- 
ficient wavelength coverage shortward of 6.0 /im, we have 
found water ice in up to 18 galaxies and upper limits of 
better than Ti CO =0.1-0.3 for the absence of water ice ab- 
sorption in another 28 galaxies. Although a small sample, 
it is of interest to correlate galaxy types with the mid- 
IR absorption/emission characteristics. The results of our 
analysis are presented in Table ||. 

Table U suggests that water ice might be a common 
species in ULIRGs, since 12 out of 19 do show ice. Seyferts, 
on the other hand, seem to be ice-poor, with ice detected 
in only 2 out of 62 galaxies. For starburst galaxies the 
numbers are less obvious, with 4 detections on a total of 
21 galaxies. 

In order to investigate this issue further we have com- 
posed average spectra for the three main galaxy types 
listed in Table |. Since the presence of 5-6 /im PAHs in 
spectra complicates the detection of the 6.0 /im water ice 
feature, we have decided to generate separate average 
spectra for continuum-dominated Seyferts and ULIRGs 
and for PAH-dominated Seyferts and ULIRGs. We clas- 
sify a spectrum as continuum-dominated when the 6.2 /im 
PAH line-to-continuum ratio is less than 0.15 and PAH- 
dominated when the ratio exceeds this value. The resulting 
five average spectra are shown in Fig. [l^ and are based on 
138 out of 250 galaxy spectra in our database. The galax- 
ies included in the averaging process were selected on the 
basis of their noise after scaling by their individual 6.5 /im 
flux. 

The top panel of Fig.|l3| shows the average PAH- 
dominated Seyfert galaxy, composed of 34 galaxy spectra. 
In addition to the well-known PAH emission features, the 
spectrum shows the 6.99 /im [Ar 11] and 10.54 /im [S iv] for- 
bidden lines, the latter of which usually is strong in active 
galaxies. The spectrum bears no trace of 6.0 /im water ice 
absorption. The second panel of Fig. |l3| depicts the aver- 
age continuum-dominated Seyfert, composed of 45 galaxy 
spectra. The spectrum is dominated by the [S iv] line, with 
traces of [Ar 11] , maybe of 7.65 /im [Ne vi] and of weak PAH 
emission — the latter in accordance with our selection cri- 
terion of 6.2/mi-PAH L/C<0.15. The 9.7/im silicate ab- 
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sorption feature is strikingly absent, as is the 6.0 /im water 
ice absorption feature. The third panel of Fig.[l3| shows 
the average continuum-dominated ULIRG, composed of 
11 galaxy spectra, among which are all Class 1 sources ex- 
cept for NGC 4418, as well as the Class 2 source Mrk 231. 
The spectrum is dominated by a broad flux peak cen- 
tered at ~8 fim. A 6.0 /im water ice feature can be easily 
recognized, starting at ~5.6/im and reaching maximum 
depth at 5.9-6.0 /im. The spectrum is noticeably differ- 
ent from the other continuum-dominated spectrum, the 
Seyfert spectrum, in the panel above. The fourth panel 
shows the average PAH-dominatcd ULIRG, composed of 
19 galaxy spectra. Among the selected sources are several 
Class 2 and Class 3 sources. The spectrum shows strong 
resemblance to the average starburst galaxy (panel below; 
29 spectra), with two exceptions: (1) the average PAH- 
dominated ULIRG shows a clear 6.0 /im water ice feature, 
partially filled in by 6.2 /im PAH emission, the average 
starburst galaxy does not; and (2) the average starburst 
spectrum shows a hint of 5-6 /im PAH features (Fig.|f), 
where the average PAH-dominated ULIRG as well as the 
other three average spectra do not. The presence of the 
5-6 /xm spectral structure near 6.0 /im forces a higher up- 
per limit to the amount of ice in the average starburst 
spectrum than in the two average Seyfert spectra. 

The average spectra largely confirm our initial find- 
ings: ULIRGs have on average more ice than Seyfert and 



Table 2. Break-down of ice galaxies over the three main 
galaxy types in our sample of 103 ISO galaxies having 
good S /N spectra and sufficient wavelength coverage blue- 
ward of 6 /im. 



Gal. type 


PHT-S 


CAM-CVF 


total 


Seyfert 


2/55 


0/7 


2/62 


ULIRG 


9/12 


3/7 


12/19 


Starburst 


3/12 


1/9 


4/21 


Other 


0/1 


0/0 


0/1 



starburst galaxies. The absence of contaminating PAH 
emission in the continuum-dominated Seyfert spectrum 
implies a tighter upper limit for the presence of water ice 
absorption than in the PAH-dominated Seyfert and star- 
burst spectra. 

Also of interest is the presence of a broad flux peak 
centered at ~8/im in the continuum-dominated ULIRG 



(Fig. 13, middle panel). In the absence of strong 6.2/im 
PAH emission (L/C<0.15), the 8/im flux peak cannot be 
attributed to strong and exceptionally broad PAH emis- 
sion. We will come back to this issue later this Section. 

The galaxies in our sample have turbulent and hos- 
tile environments (i.e., the diffuse interstellar medium or 
the immediate surroundings of the central starburst or 
AGN) which significantly process and often destroy the 
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Fig. 13. Panels showing the average spectra of PAH- 
dominated (6.2^m-PAH L/C>0.15) and continuum dom- 
inated (6.2/Ltm-PAH L/C<0.15) Seyferts, ULIRGs and 
starburst galaxies. Grey shaded areas indicate the stan- 
dard deviations in the average spectra. Vertical dotted 
lines denote the central wavelengths of well-known spec- 
tral features (see text). 



molecular material that reside in or close to these re- 
gions. However, within our ice-galaxy sample, volatile ma- 
terial (which is easily destroyed) is clearly present. This 
suggests that some shielding from harsh environments 
must occur. In the diffuse medium supernova shock waves 
are the predominant destruction mechanism for interstel- 
lar dust (Jones et al. 



1994) 



Refractory grains and most 
likely PAHs are typically destroyed within 10 s yr, whereas 
ice grains of 1000 A can only survive for ~10 6 yr. These 
lifetimes increase significantly in dense clouds where the 
molecules are protected from processes such as sputter- 
ing by shocks. On the other hand, molecular material in 
an AGN torus (or in nearby dusty clouds) is exposed to 
the intense hard X-rays originating from the central en- 
gine. Due to their extreme energies, hard X-ray photons 
are able to penetrate column densities of ~ 10 25 cm~ 2 (e.g. 



NGC4945: Iwasawa et al. 1993 ). Hence, under such con- 
ditions, ice mantles on the dust grains would be unable 
to survive. PAHs, on the other hand, might be able to 
survive for lifetimes of > 10 4 yr if the column density is 
greater than ~ 10 23 cm~ 2 (Voit 1992) and hence it is more 



likely to see PAHs as opposed to ices in such environ- 
ments. Nevertheless, the harsh conditions in AGN toroids 
may well influence the composition of the emitting PAH 
family. Further modelling and laboratory experiments are 
required to determine to what extent PAHs and ices are 
processed by the hard radiation present in AGNs. 

ULIRGs are merger systems in which the tidal forces 
have allowed molecular material from the galaxy disks to 
accumulate in the nuclei of the system (e.g. Solomon et 
al. 119971; Downes & Solomon tL998t Sakamoto et al. 11999; 



Tacconi et al. 1999). The presence of water ice in ULIRGs 



thus comes as no surprise. In the many dense shielded 
cloud cores, that are likely to be present, water ice can 
survive until these clouds are dispersed by starformation. 
Clearly, the "icy" ULIRGs still have ample cold dense 
molecular material to continue star formation for some 
time. In Galactic dense molecular clouds the 6.0 fj.ni wa- 
ter ice absorption feature is seen against the mid-IR con- 
tinuum of deeply embedded protostars, like W33A and 
NGC 7538 IRS9. In external galaxies, the background con- 
tinuum source might also be provided by hot dust associ- 
ated with the AGN torus or with starburst activity. The 
similarity of the average ULIRG spectrum (Fig. 1 of Lutz 
et al. 1998) to Class 3 spectra with strong PAHs, suggest 
that ice absorption in ULIRGs is preferentially linked to 
a scaled up star formation process. The lower luminos- 
ity example of NGC 4418 cautions, however, that this link 
cannot be generalized to all individual objects: ice absorp- 
tion may also be seen against AGN. 

Little can be derived from our data on the ice ab- 
sorption properties of lower luminosity starbursts. The 
very strong PAH emission features prevent the detection 
of moderately strong 6.0 /im features. In fact, given the 
lower obscuration of starbursts compared to ULIRGs or 
to Galactic protostars, weaker ice features in starbursts 
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Table 3. Overview of features observed in Class 4-6 
galaxies shown in Figs.|[ ||&[l(], all displaying no sign 
for 6/im water ice absorption at the r(6/imice)=0.1-0.3 
level. 



Target 


Galaxy 


cz 


5 /im 


9.7 /im 




type 




PAH 


silicate 






[km/s] 


emission 


absorption 


NGC 253 


SB 


245 


++ 


? 


NGC 520 


SB 


2281 


++ 


+ 


NGC 613 


SB 


1475 




? 


NGC 1068 


Sy2 


1148 




+ 


NGC 1365 


Syl.8 


1636 


? 


? 


104385-0828 


Sy2 


4527 




++ 


NGC 1808 


Sy2 


1000 


++ 


? 


MGC 8-11-11 


Syl.5 


6141 






M82 


SB 


203 


++ 


? 


M83 


SB 


516 


++ 


? 


Mrk279 


Syl.5 


8814 






Circinus 


Sy2 


439 


? 


++ 


PKS 2048-57 


Sy2 


3402 




+ 



++ = strongly present 
+ = present 
? = maybe present 
- = absent 



are expected for similar properties of the obscuring mat- 
ter. The examples of M 82 and NGC 253 (Sturm et al. 
200C ) suggest that the 3 /im region may be better suited 
to quantify such weaker ice absorptions, with the caveat 
for ground-based observations of the 3 /im region being a 
difficult atmospheric window. Ice is seen in these starbusts 
at 3 /im but not detectable in the heavily structured 6 /im 
spectrum. 

The conditions for the existence of water ice near 
Seyfert nuclei seem to be generally poor. If most of the ob- 
scuration occurs in a dense torus that is surrounding the 
nuclear region and heated by the X-rays from the central 
engine, ices are unlikely to survive. Indeed, 3 /im spectra 



of AG N (Bridger et al. |1994| ; Wright et al. |199J; Imanishi 
2000| ; Sturm et al. |200C| ) tend to show 3.4 /im hydrocarbon 
absorptions but no ices. 

The low number of Seyfert and starburst ice-galaxies, 
as inferred from the mid-IR spectra, might indicate that 
dense shielded molecular cloud cores are far less abun- 
dant in these systems than in ULIRGs. Likely, most of 
the absorption in Seyferts originates in the dense toroid 
surrounding the nuclear region, which is thoroughly en- 
ergetically or thermally processed by the X-rays of the 
internal monster and therefore contains no ices. 

As a corollary and as illustrated by Fig.|^ the presence 
of silicates does not automatically imply the presence of 
water ice. If the environment is not sufficiently shielded or 
is too warm, water ice cannot exist. For the five galaxies 
displayed (one starburst galaxy and four Seyferts) this 
might be the case. 



Conversely, does water ice correlate with the presence 
of shielding dust (e.g., silicates) ? Except for 105189-2524, 
all our ice-galaxies show a strong 9.7 /xm silicate absorp- 
tion feature (see Table]!]), with t s h >1.5 for most sources. 
The ratio Ti ce /r S ii spans a range of <0.10 for 100183-7111 
to 1.8 for 123128-5919, with a mean value of 0.67 for the 
galaxies with a non-saturated silicate feature. In contrast, 
all five Class 1 galaxies have Ti cc /T s u <0.37. This suggests 
that silicate and water ice optical depths are not corre- 
lated on the galaxy-wide scale that we probe with our 
ISO observations. We note that low mass protostars in in- 
dividual clouds (e.g. Taurus) do show a good correlation 
between ice optical de pth and the strength of the silicate 
feature (Whittet et al. |l983| , |l988[ ). However, for massive 
luminous protostars in different clouds, no such correlation 
exists, largely due to the extensive thermal processing of 
the environment by the newly formed star (cf. Tielens & 
Whittet 1997 ). Obviously, the abundance of water ice is 
controlled by special conditions, which in our sample are 
apparently only met in some of the regions within our 
Class 1 ice galaxies. 

Though originally published as a unique source, it 
is apparent from our study that some characteristics of 
NGC 4418 are seen in a larger number of galaxies. Mid-IR 
absorptions of water ice are detected in about one sixth 
of the galaxies for which ISO spectra of adequate S/N 
are available, and perhaps the majority of ULIRGs. These 
detections range from spectra dominated by an absorbed 
continuum (Class 1) over spectra mainly showing absorp- 
tion signatures but also clear PAH emission (Class 2) to 
spectra that resemble normal PAH emission spectra sub- 
jected to additional dust and ice obscuration (Class 3). 

Short of full spectral modelling, empirical criteria help 
to assess the relative importance of emission and absorp- 
tion in shaping these spectra. First and most important 
is the distinction between a 6.2 /im PAH emission and 
the 6.5-6.7 /tm pseudo maximum in an absorption dom- 



inated spectrum (Spoon et al. 2001). Second, there are 
differences in the shape of the maximum in the 8 /im re- 
gion. PAH dominated spectra exhibit rel ativel y narrow 7.7 
and 8. 6 /im emissions (e.g. Helou et al. 2000 Rigopoulou 



et al. 1999|) . Also, they show a very steep rise between 



7.2 and 7.7/tm. Obscuration of such a spectrum (assumed 
to be intrinsically unchanged) affects mainly the relative 
strengths of such features but keeps the 7.7 /im peak nar- 
row («0.6/im FHWM, e.g. Fig.0; Fig. 6 of Rigopoulou et 



al. 1999 ) and its ascent steep. In contrast, absorption dom- 
inated spectra tend to show a wider, less well defined peak 
around 8 /tm (e.g. Fig.|i~I|), with significant variation de- 
pending on continuum shape and optical depth of the var- 
ious features. 

The change from an absorbed continuum to an ab- 
sorbed PAH spectrum is illustrated further in spectra 
of three adjacent positions in the Galactic star forming 
complex W3 (Fig. [l^ , D. Cesarsky, priv. comm.). The 
top panel towards the infrared source is dominated by 
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a continuum absorbed by ice, HAC, and silicates with a 
broadish 8 /im bump. The middle panel is an intermediate 
case and the bottom shows an obscured PAH spectrum. 
These spectra illustrate that the qualitative features dis- 
cussed above and outlined in the schematic views of Fig. || 
and Fig.[l4| can indeed originate from radiation transfer 
through a dense and dusty medium, and that star forma- 
tion in our own galaxy can reproduce the full range of 
phenomena. 

Differences in the shape of the 8 /im maximum may 
help in understanding the vast majority of sources that 
are not pure absorbed continua. Is the frequently observed 
weakness of the 6.2 /im feature (compared to the 7.7/8 /im 
maximum) due to obscuration of a PAH-dominated spec- 
trum or due to superposition of a little obscured PAH 
spectrum and a strongly obscured continuum? The intrin- 
sic 6.2/7.7 PAH ratio is qu ite stable for many Galactic 
sources (Peeters et al. 2002| ) as well as the ISM in normal 
and starburst galaxies (Helou et al. 200C; Rigopoulou et 
al. |1999|, ou r Fig.|). Lutz et al. ( [1998; ) and Rigopoulou 
et al. (1999) argued in favour of obscuration as cause of 
the 6.2 weakness, using an observed correlation between 
6.2/7.7 feature ratio and extinction to the starburst re- 
gion, as measured from independent mid-IR emission line 
data (Genzel et al. 199S). On the other hand, suggestions 
have been made of an effectively unobscured 'surface layer' 
producing most of the PAH emission, CII line emission, 



and submm continuum of ULIRGs (e.g. Fischer et al. 2001 
Haas et al. ^001| ). Such an approach tackles problems like 



the observed [C n] deficit of ULIRGs at the expense of 
breaking the link between PAHs and the star formation 
observed in the mid-infrared fine structure lines, and of 
introducing an unknown fully obscured component. 

Water feature optical depths of r ico ~1 (Table)!]) are 
sufficient to modify the ratio of the 6.2 and 7.7 /im PAH 
features by a factor ~2. While this is enough to explain 
the weakness of the 6.2 /im feature in many objects by ob- 
scuration, both the uncertainty in the optical depth values 
and the lack of knowledge to which extent the ice features 
apply to the PAHs and/or an underlying continuum make 
conclusions for individual objects very difficult. Here, use 
of the shape of the feature near 7.7/8 /im can help with 
breaking the degeneracy. As long as the intrinsic shapes 
of the PAH features are assumed to remain constant, a 
wider than usual peak near 8 /im argues for a consider- 
able contribution of a heavily absorbed continuum. Our 
Class 1 and 2 spectra (Fig.[| and ^) suggest that this is 
the case for a number of galaxies. In a future paper (Spoon 
et al. in prep.) we will derive quantitative fits to the high 
quality Class 2 spectrum of Arp 220, arguing for a best fit 
with a considerable contribution by an absorbed contin- 
uum. The existence of Class 1 and Class 2 sources implies 
that estimates of the PAH contribution on the basis of a 



ably for the simple line-to-continuum ratio. This effect is 
highly relevant for some of the Class 1 and 2 spectra, but 
less so for Class 3 or the ULIRGs in general. High S/N 
spectra, extended wavelength coverage, and fits of the en- 
tire wavelength range with PAHs and continua obscured 
by dust and ices are needed for quantitative progress. 

The issue of broad 8 /im features is complicated further 
by the presence of yet another category of " 8 /im maxima" 
in addition to PAHs and absorbed spectra: Several lumi- 
no us AG N-like ULIRGs (Tran et al. |200l| ; Taniguchi et 
al. 1997) show broad 8 /im maxima on top of a smooth 



mid-infrared continuum, with PAH absent or very weak. 
Clear examples for this include 109104+4109, 100275- 
2859, 122192-3211, and 123529-2119. The key difference 
to our Class 1 and Class 2 spectra with broad 8 /im max- 
ima is that these objects do not show the very deep silicate 
feature of the absorbed Class 1 and 2 spectra. Tran et al. 



(2001) discussed these features in terms of self-absorbed 
silicate emission or a modified PAH origin, without defi- 
nite c onclus ion. Recent survey s of g alactic PAHs (Peeters 
et al. 2002: Vcrstractc ct al. |200l|) fail to observe simi- 



lar profiles even under unusual conditions. The nature of 
these features remains uncertain. 

It is interesting to speculate in analogy to Galactic 
sources that our classification might reflect an evolution- 
ary sequence. If an evolved starburst is represented by a 
a PAH spectrum, while strong mid-infrared continua are 
typical for the H H-regio n con tinuum of intense compact 
starbursts (Laurent et al. 2000 ) or AGN, then Class 1 may 
trace the deeply embedded beginnings of star formation 
(or AGN activity), while the latter classes reflect more 
advanced and less enshrouded stages. 

Finally, a point of caution. Since the 3.0 /tm water ice 
absorption extends over a difficult part of the L band, its 
definition is not easy in ground based spectra like those of 
Imanishi et al. (2001). Its short end is outside the L band 
and its long end coincides with the 3.3 /im PAH emis- 
sion and 3.4 /im hydrocarbon absorption. A gradually in- 
creasing continuum from 3.15 to 4.0 /im is however visible 
in their spectrum of UGC5101. Depending on the loca- 
tion along the line of sight, the water ice absorption will 
strongly affect measured quantities like the 3.3 /im contin- 
uum, the 3.3 /im PAH flux, its equivalent width and the 
3.3/im-PAH-to-far-IR flux ratio. Hence, in galaxies with 
water ice absorption (like UGC5101, Tj cc =1.3), both un- 
corrected 3.3 /im and 6.2 /im PAH fluxes will provide un- 
derestimates when used as star formation indicators. For 
the 3.3 /im PAH, aperture losses of present instruments 
are an additional reason for underestimates as in the case 
of UGC5101 where the north-south oriented 1.2" slit of 



Imanishi et al. (2001) is perpendicular to the orientation 



of the >2" star forming region (Genzel et al. 1998) 



line-to-continuum ratio (Genzel et al. 1998) or of fit s with 



extinction laws that do not include ices (Tran et al. 2001 ) 
will overestimate the importance of PAHs, more notice- 
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7 8 9 

Rest wavelength [/im] 
Fig. 14. Schematic view of the impact of dust and ice 
absorption on a mid-IR continuum spectrum. All shaded 
areas combined constitute the adopted local mid-IR con- 
tinuum. The black area is all that is left of the local con- 
tinuum after ices and silicates absorb substantial portions 
of the 5-12 /im local continuum. Note the presence of an 
emission-like broad feature near 8 /im in the resulting ob- 
served spectrum, which at first glance or at poor S /N may 
be mistaken for 7.7 /im PAH emission. 



7. Conclusions 

Following the discovery of water ice in NGC 4418 (Spoon 
et al. EOOlD , we have searched our sample of ISO galaxy 
spectra for galaxies showing similar signs of 6.0 /im water 
ice absorption. We have found in total 18 galaxies, which 
we grouped into three classes. The classification is based 
on the presence of ice absorption, 6.2 /im PAH emission, 
and the nature of the 7.7-8 /im feature: PAH emission or 
absorbed mid-IR continuum, or a combination of the two. 
We also looked for galaxies showing no signs of water ice 
absorption. This sample contains 28 galaxies with water 
ice upper limits ranging from rj co =0.1-0.3, depending on 
the S /N and the complexity of the 5-7 /im spectrum. We 
classified these galaxies in another three classes (Class 4- 
6), depending on the presence of 9.7 /im silicate absorption 
and 5-6 /im PAH emission. 

We have modeled the complicated interplay of 6.0 /im 
water ice absorption and 6.2 /im PAH emission, which in 
ISO-PHT-S spectra takes place at the blue end of the 
PHT-SL range and, hence, only can be recognized in 
galaxies with a redshift in excess of 3000 km/s. For spec- 
tra obtained with the ISO-CAM-CVF instrument no such 
limitations apply. Our modeling supports the presence of 
ice in the Class 1-3 sources, except in a few cases where 
ice is not strictly required to obtain an acceptable fit to 
the observed spectrum. 

Based on a subsample of 103 good S/N ISO galaxy 
spectra with sufficient wavelength coverage blueward of 




7 8 9 10 11 
Rest wavelength [/im] 

Fig. 15. ISO-CAM-CVF spectra of three adjacent po- 
sitions in the W3 molecular cloud (D. Cesarsky, priv. 
comm.), showing a gradual transition between an ob- 
scured PAH-dominated mid-IR spectrum (bottom) and 
a strongly absorbed hot dust continuum (top). The ver- 
tical dotted lines show the positions of 6.0 /mi water ice, 
6.25 /im PAH, 6.85 /im HAC and 7.7 /im & 11.3 /im PAH. 



6 /im, a substantial fraction of ULIRGs (12 out of 19) 
contain detectable amounts of water ice. On the other 
hand, the majority of Seyfert (2 out of 62) and starbursts 
galaxies (4 out of 21) probably do not. These results are 
confirmed by the spectral structure seen in our average 
spectra of Seyfert, ULIRG and starburst galaxies: water 
ice absorption is only obvious in the average ULIRG. 

Class 1 & 2 ice galaxies are dominated by a broad fea- 
ture near 8 /im which indicates a strong contribution by 
a dust and ice-absorbed continuum, similar to that seen 
in NGC 4418. These observation stress the need for high 
S/N data and refined diagnostic methods, to properly dis- 
criminate spectra dominated by PAH emission and spectra 
dominated by heavy obscuration. 

The interplay of the broad 8 /im feature and PAH emis- 
sion, as seen in our ice galaxies, shows strong similarities 
with features seen in Galactic star forming clouds. This 
leads us to believe that our classification of ice galaxy spec- 
tra in three classes might reflect an evolutionary sequence 
from strongly obscured beginnings of star formation (and 
AGN activity) to a less enshrouded stage of advanced star 
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formation (and AGN activity), as the PAHs get stronger 
and the broad 8 /im feature weakens. 

Acknowledgements. The authors wish to thank Diego Ce- 
sarsky, Els Peeters and Eckhard Sturm for discussions. This 
research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Lab- 
oratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 

References 

Allamandola L.J., Bregman J.D., Sandford S.A., et al., 1989, 
ApJ 345, L59 

Bohlin R.C., Savage B.D. & Drake J.F., 1978, ApJ 224, 132 
Bridger A., Wright G.S. & Geballe T.R., 1994, in Infrared As- 
tronomy with Arrays, I. McLean (ed.), Kluwer-Dordrecht, 
537 

Cesarsky C.J., Abergel A., Agnese P., et al., 1996, A&A 315, 
L32 

Chiar J.E., Tielens A.G.G.M., Whittet D.C.B., et al., 2000, 
ApJ 537, 749 

Clavel J., Schulz B., Altieri B., et al., 2000, A&A 357, 839 
Dartois E. & d'Hendecourt L., 2001, A&A 365, 144 
De Graauw Th., Haser L.N., Beintema D.A., et al., 1996, A&A 
315, L49 

Downes D. & Solomon P.M., 1998, ApJ 507, 615 

Draine B.T., 1989, in Infrared Spectroscopy in Astronomy, B.H. 

Kaldeich (ed.), ESA-SP 290, 93 
Fischer J., 2001, in The promise of the Herschel Space Observa- 
tory, ESA-SP460, G.L.Pilbratt, J. Cernicharo, A.M. Heras, 
T. Prusti, R.Harris (eds.), ESA-SP460, 131 
Furton D.G., Laiho J.W. & Witt A.N., 1999, ApJ 526, 752 
Genzel R., Lutz D., Sturm E., et al., 1998, ApJ 498, 579 
Haas M., Klaas U., Miiller S.A.H., Chini R. & Coulson I., 2001, 
A&A 367, L9 

Hagen W., Tielens A.G.G.M. & Greenberg J.M., 1983, A&AS 
51, 389 

Helou G., Lu N.Y., Werner M.W., Malhotra S., et al., 2000, 
ApJ 532, L21 

Hudgins D.M., Sandford S.A., Allamandola L.J. & Tielens 

A.G.G.M., 1993, ApJS 86, 713 
Imanishi M., 2000, MNRAS 319, 331 

Imanishi M., Dudley C.C. & Maloney PR., 2001, ApJ 558, L93 
Iwasawa K., Koyama K., Awaki H., et al., 1993, ApJ 409, 155 
Jones A., Tielens A.G.G.M., Hollenbach D.J. & McKee C.F., 

1994, ApJ 433, 797 
Keane J.V., Tielens A.G.G.M., Boogcrt A.C.A., et al., 2001, 

A&A in press 

Kessler M.F., Steinz J.A., Anderegg M.E., et al., 1996, A&A 
315, L27 

Laurent O., Mirabel I.F., Charmandaris V., et al., 2000, A&A 
359, 887 

Lemke D., Klaas U., Abolins J., et al., 1996, A&A 315, L64 
Lutz D., Spoon H.W.W. , Rigopoulou D., et al., 1998, ApJ 505, 
L103 

Lutz D., 1999, In: The Universe as seen by ISO, P. Cox, M.F. 

Kessler (eds.), ESA-SP 427, 623 
Maldoni M.M., Smith R.G., Robinson G. & Rookyard V.L., 

1998, MNRAS 298, 251 
Moorwood A.F.M., Lutz D., Oliva E., et al., 1996, A&A 315, 

L109 



Moutou C, Sellgrcn K., Leger A., et al., 1998, in Star For- 
mation with the Infrared Space Observatory, J.L. Yun, R. 
Liseau (eds.), ASP Conf. Series 132, 47 

Peeters E., Hony S., Van Kerckhoven C, et al., 2002, A&A 
accepted 

Pendleton Y.J., Sandford S.A., Allamandola L.J. & Tielens 
A.G.G.M., 1994, ApJ 437, 683 

Pendleton Y.J. & Chiar J.E., 1997, in The Nature and Evo- 
lution of Interstellar Organics From Stardust to Planetes- 
imals., Y.J. Pendleton & A.G.G.M. Tielens (eds.), ASP 
Conf. Ser. 122, 179 

Rigopoulou D., Spoon H.W.W., Genzel R., et al., 1999, AJ 
118, 2625 

Roche P.F. & Aitken D.K., 1984, MNRAS 208, 481 
Sakamoto K., Scoville N.Z., Yun M.S., Crosas M., et al., 1999, 
ApJ 514, 68 

Solomon P.M., Downes D., Radford S.J.E. & Barrett J.W., 

1997, ApJ 478, 144 
Spoon H.W.W., Koornneef J., Moorwood A.F.M., et al., 2000, 

A&A 357, 898 

Spoon H.W.W., Keane J.V., Tielens A.G.G.M., Lutz D. & 

Moorwood A.F.M., 2001, A&A 365, L353 
Sturm E., Lutz D., Tran D., et al., 2000, A&A 358, 481 
Tacconi L.J., Genzel R., Tecza M., Gallimore J.F., Downes D., 

Scoville N.Z., 1999, ApJ 524, 732 
Taniguchi Y., Sato Y., Kawara K., et al., 1997, A&A 318, LI 
Tielens A.G.G.M. & Whittet D.C.B., 1997, in Molecules in As- 
trophysics: Probes and Processes, E.F. van Dishoeck (ed.), 
Kluwer-Dordrecht, 45 
Tran Q.D., Lutz D., Genzel R., et al., 2001, ApJ 552, 527 
Verstraete L., Pech C, Moutou C. et al., 2001, A&A 372, 981 
Voit G.M., 1992, MNRAS 258, 841 

Whittet D.C.B., Bode M.F., Baines D.W.T., et al., 1983, Na- 
ture 303, 218 

Whittet D.C.B., Bode M.F., Longmore A. J., et al., 1988, MN- 
RAS 233, 321 

Whittet D.C.B., Schutte W.A., Tielens A.G.G.M., et al., 1996, 
A&A 315, L357 

Wright G.S., Bridger A., Geballe T.R. & Pendelton Y.J., 
1996, in Cold Dust and Galaxy Morphology, D.L. Block, 
J. M. Greenberg (eds.), 143, Kluwer-Dordrecht 



H.W.W. Spoon et al.: Ice features in the mid-IR spectra of galactic nuclei 




Fig. 8. Class 3 (continued) 
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Fig. 9. Class 4 mid-IR galaxy spectra. Galaxies of this kind are ice-poor and show the presence of silicates along the 
line of sight. The presence of silicates is inferred from the absorption feature at 9.7 /im, characteristic for silicates. The 
upper limit on H 2 ice is inferred from the absence of a 6.0 /im absorption feature. For NGC 1068 and for Circinus 
three spectra are shown: the ISO-SWS nuclear spectrum in grey and the ISO-CAM-CVF nuclear and total spectra in 
black. The offset between the three spectra is most likely an aperture effect. The panels on the left show the 5-12 ^m 
spectra on a linear scale, while the panels on the right show the same spectra on logarithmic scale, with smoothed and 
rcbinncd 2-5/im ISO-PHT-S spectra and 12-100 /im IRAS fluxes added. All fluxes have been k-corrected. The vertical 
error bars at 5.8 and 11 /im denote the pixel to pixel flux error (la, so not ±lc), while the horizontal error bars denote 
the gradually changing ISO-CAM-CVF spectral resolution (R~40). 
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Fig. 10. Class 5 mid-IR galaxy spectra. Galaxies of this kind are ice-poor and show no clear evidence for the presence 
of silicates along the line of sight. Note the clear mismatch of the IRAS FSC fluxes with the ISO-CAM-CVF nuclear 
spectra for both NGC613 and NGC 1365. This most probably is an aperture-size effect. The presence of strong PAH 
emission both short- and longward of 9.7 /mi makes it difficult to completely rule out silicate absorption in the nuclei 
of NGC 613 and NGC 1365. Compared to Circinus and NGC 520 however (Fig.^), the silicate depth must at least be 
small. The panels on the left show the 5-12 /tm spectra on a linear scale, while the panels on the right show the same 
spectra on logarithmic scale, with smoothed and rebinned 2-5/im ISO-PHT-S spectra and 12-100 /im IRAS fluxes 
added. All fluxes have been k-corrected. The vertical error bars at 5.8 and 11 /tm denote the pixel to pixel flux error 
(la, so not ±1<t), while the horizontal error bars denote the gradually changing ISO-CAM-CVF spectral resolution 
(R~40). 
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Fig. 11. A comparison of all the sources with the best fitting models which contain both H2O ice absorption and PAH 
emission bands. The grey histogram data represent the model fits. The three sources plotted in Fig.Jl^ are not shown 
here. 
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Fig. 12. Examples of 3 sources which show only a few pixels in the 6.0 fj,m range and which might hinder the 
interpretation of H2O ice being present. The left panel shows the best fits by a model containing both H2O absorption 
and 5-6 /jm PAH emission. The right panel shows the best possible fits by a model which only contains 5-6 /im PAH 
emission features. The adopted continua are denoted by the dashed lines and the resulting fits are represented by 
the line shown in grey. Note that since our model does not take into account the 7.7 and 8.6 /mi PAH features (see 
Sect. 3.1), the quality of the fit can only be tested shortward of 7/xm. 



